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Abstract 

Here we present the program for the transformation of the expansions of the atomic state functions 
from LS- to jj-coupled bases. The program is a part of the coupling optimization package Laris and 
allows to transform the multiconfigurational expansions with the arbitrary number of open shells (includ- 
ing /-electrons) and to identify atomic levels and classify the atomic spectra by means of intermediate 
quantum numbers of LS- and jj-couplings. 
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PROGRAM SUMMARY 



Title of program: LSJJ 

Program obtainable from: Vilnius University Research Institute of Theoretical Physics and Astronomy, A. Gostauto 
12, Vilnius, 2600, Lithuania. E-mail: tomas@mserv.itpa.lt 

Computer for which the library is designed and others on which it has been tested: Pentium-based PC 155 MHz. 
Installations: Vilnius University Research Institute of Theoretical Physics and Astronomy (Lithuania) 
Operating systems or monitors under which the new version has been tested: LINUX 2.4.5 

Programming language used in the new version: ANSI standard Fortran 90/95. 

Memory required to execute with typical data: Memory requirements depend on the shell structure and the size of the 
wave function expansion which is used to represent the atomic levels. 

No. of bits in a word: All real variables are parametrized by a selected kind parameter and, thus, can easily be 
adapted to any required precision as supported by the compiler. Presently, the kind parameter is set to double 
precision (two 32-bit words) in the module rabs_constant. 

No. of bytes in distributed program, including test data, etc.: ... 

Distribution format: compressed tar file 

Keywords : atomic structure, LS*-coupling, jj-coupling, LS — jj transformation, optimal coupling scheme, classifica- 
tion of atomic spectra. 

Nature of physical problem: 

The set of quantum numbers for the identification of atomic state functions (levels) and the classification of atomic 
spectra usually are chosen comparing the distributions of the mixing coefficients in the bases corresponding to various 
coupling schemes. The corresponding mixing coefficients usually are obtained diagonalizing the energy matrix in the 
basis of the wave functions of certain coupling. Once the secular equation is solved and the mixing coefficients in 
some coupling scheme are obtained they can be used as a source to get corresponding coefficients for the bases of 
other coupling schemes if the appropriate transformation matrix is known [1] . 

Method of solution: 

For the transformation of the atomic state functions from LS- to jj-coupled basis we use the recently obtained LS~jj 
transformation matrices for the shells of equivalent electrons [2] as well as the formula for the LS — jj transformation 
matrices [3]. 

Restrictions on the complexity of the problem: 

The program restricted to the configurations with LS-coupled shells with 1=0, 1, 2 and 3. 

Unusual features of the program: The LSJJ program is designed as a part of the coupling optimization package Laris. 
References: 

[1] Z. Rudzikas, Theoretical Atomic Spectroscopy (Cambridge Univ. Press, Cambridge, 1997). 
[2] G. Gaigalas, T. Zalandauskas and Z. Rudzikas, At. Data Nucl. Data Tables, in print (2002). 

[3] G. Gaigalas, Z. Rudzikas and T. Zalandauskas, Book of Abstracts of XI HCI Conference, 1-6 September, 2002, 
Caen, France, 48. 
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LONG WRITE-UP 



1 Introduction 

In investigating the energy spectra and other characteristics of atoms and ions, the experimentalists most 
often use the LS'-coupling scheme. There a shell of electrons is made of the electron states having the same 
principal and angular quantum numbers, n and I. Such a shell is additionally characterized by the shell 
occupation number N and the resulting orbital and spin momenta, Li and Si. When several shells are 
present, the state in LS-coupling is additionally characterized by the intermediate resulting momenta L^ 
and Sij, and by the total momenta L, S and J. This coupling scheme is very popular in the theoretical 
atomic physics when atoms and ions are investigated in the non-relativistic approximation. This, in the first 
place, is related to the fact that for the light atoms and ions such an approximation produces quite accurate 
results and the agreement with the experiment is fairly good. Secondly, this coupling scheme is realistic for 
such systems. Therefore it is convenient for the identifying the energy levels. At the same time, this scheme 
is convenient to investigate such processes as transitions, whereby the selection rules for electronic transitions 
are both easily determined and reflect the real situation. But this coupling scheme is for a relatively small 
number of atoms and ions. The quantum numbers of this scheme are not accurate for highly charged ions 
and a number of other cases. For the identification of atomic state functions (ASF) and classification of 
energy spectra the approximate quantum numbers are widely used. Usually such a classification relies on 
the characteristics of the LS— or jj— coupling scheme. But in a number of cases (for example ions of 
intermediate ionization degree) the choice of such classification is not trivial and implies the consideration 
of the expansions of atomic state functions in several bases of pure coupling wave functions. Therefore it 
is convinient to have the possibility to transform the expansions of atomic state functions between various 
bases of pure coupling wave functions [1]. 

The transformation matrices between the different bases of pure coupling wave functions were studied 
([2]-[10]), whereas the corresponding software were developed ([11, 12, 13, 14]). Nevertheless until recently 
there were neither LS — jj transformation matrices for the arbitrary number of open shells available, nor 
the corresponding software allowing one to transform the ASFs resulted from the large scale calculations. 
In addition, the software created earlier does not always allow to find the most suitable coupling scheme, 
as sometimes it appears that no dominating scheme is present and several schemes are equally suitable 
simultaneously. We present the investigation of the Fluorine isoelectronic sequence as an example [15]. In 
this case, several coupling schemes are treated in the search for an optimal one. Those are LS, LK, JK 
and J J. In all those schemes the LS coupling is valid inside a shell. Using the analysis of the structure of 
the weights of the ASF or the quantity R, one may notice that in this particular isoelectronic sequence the 
ions do not have a dominating scheme, or the dominating scheme does not stand out clearly among others. 
Therefore the above-mentioned methods have been extended [8, 9, 10] in a way providing a possibility to 
change the LS— coupling scheme into the jj— coupling one inside a shell, when necessary. The recently 
obtained LS — jj transformation matrices for the shell of equivalent electrons [9] as well as the formulas for 
the LS — jj transformation matrices [10] enabled us to develop a computer code allowing to transform the 
ASFs from LS- to jj-coupling for all the configurations practically needed in atomic physics. In the search 
for an optimal coupling scheme that would allow one to change the coupling scheme inside a shell and to 
investigate the atoms and ions with ASFs containing open /-shells with arbitrary occupation numbers. 

To facilitate such a classification of the atomic levels and wave functions, we present here the program 
Lsjj. It supports the performance of the LS — jj transformations. This program may be of value as a 
separate module in solving the physical problems that also need this transformation. However in this case 
one should bear in mind that the above-mentioned transformation would provide the correct results only if 
the calculations are performed in the framework of the Racah-Fano phase system [16], and the coefficients 
of fractional parentage are defined via the quasispin formalism, as described in [17, 18]. 

A short explanation of the theoretical background is presented in Section 2. The explanation of the usage 
and the organization of the program with the description of main subroutines are presented in Section 3. 
The example of the usage of the program and the concluding remarks are presented in Sections 4 and 5 
correspondingly. 
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2 Theoretical background 



While investigating theoretically various characteristics of atoms and ions (for example using multiconfigu- 

ration or configuration interaction approaches) the atomic state functions are used which are either written 
in terms of a LS— coupled 

l*-( jP )) = E fl - L5) ( T ) \lrL r S r JP) (1) 
r 

or j j-coupled basis 

l*^ p )) = E fl ' ],) w \^ Jp ) ' ( 2 ) 

s 

where t = 1, 2, ... enumerates the atomic states (of the given J and parity P) and j r , j s denote the sets of all 
additional quantum numbers of the LS- and jj-couplcd configuration state functions (CSF) correspondingly. 

When a pure LS'-coupling exists in the atoms or ions under consideration, then r=l in equation (1). 
When the jj-coupling exists, then s—1 in equation (2). Therefore such an ASF is identified by an array 
of quantum numbers 71. However, these quantum numbers are exact only for the cases of pure coupling 
schemes, which is more the exception than a rule. 

In the most general case the exact ASF is made up as an infinite sum in (1) or (2). Therefore an ASF 
consists of an infinite number of CSFs. In the real calculation one attempts to include terms as many as 
possible in the sum having the maximum coefficients ar LS \r) or aj^ (r). But in this case it is not clear what 
array 7 to ascribe in identifying the ASF. Bearing in mind that the coefficient modulus squared (weight) 

ar LS \r) or aj 5 ^ (r) defines the probability for the atomic state to be in a state with the array of 
quantum numbers ~f r L r S r or j s , one can easily classify all the energy levels by ascribing the characteristics 
of the maximum coefficient (ar^fi(r) in LS or ai^xC 7 ") m Ji- cou P un g) to the whole ASF. 

The situation becomes more complicated when the expansion (1) or (2) has no obviously dominant weight 
but several approximatelly equal ones. Most efficient way for such a classification is the transformation of 
the ASF into another coupling scheme with the suitable distribution of weights (i.e. with only one dominant 
weight coefficient). For example, the classification of the spectra of complex atoms and ions using the 
intermediate quantum numbers of LS-coupling quite offen is problematical. 

One of the possible solutions is to transform the function (1) into function (2). Then the level is ascribed to 

the characteristics 7r max that have the maximum weight (r) at the transformation of CSFs |7r max JP), 

i.e. the intermediate quantum numbers of jj-coupling are used for the identification of the ASFs. Such 
transformations and identifications are the tasks of the LSJJ program presented in this paper. 

For practical purposes and especially for an efficient transformation the program treats the CSF in a 
standard order, i.e. it is assumed that the LS'-coupled functions are couped consequently 

| 7r LS JP) = [(...(((/^aiLiSi, l 2 2 a 2 L 2 S 2 )L 12 S 12 , l^a 3 L 3 S 3 )L 123 S 123 )...)LS j) (3) 
as well as the jj-coupled CSFs 

hsJP) 

Ji_ - + N 1+ + _N 2 _ - + N 2+ + _N 3 _ - + N 3+ + \ 

(...(((((«! VlJl, «l flJl)Jl, K 2 V 2 J 2 )J l2 , K 2 IS 2 J 2 )J 12 K 3 ^3^3)^123, «3 V 3 J 3 )Jl 23 ...)J j (4) 

where k is the relativistic (angular momentum) quantum number Ki— k and k^= — Z, — 1. 

If both subshells with common U, i.e. Ki and Ki appear in the expansion, these two subshells always occur 
successively in the sequence (« i Vi J i7 K i v-iJi) Ji. Formally, we can use this sequence even for the subshcll 
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states with zero occupation if we interpret \k°u = J = 0) = 1; in this case, the full Clebsch-Gordan 
expansion remains valid due to the orthonormality properties of the Clebsch-Gordan coefficients. 

For the LS — jj transformation of the configuration states we further assume in a standard order that 

li =ki= — (/ti +1), I2 = k-2 = — («2 +1), i.e. that the sequence of (sub-)shell states is the same on 
both sides of the transformation matrix. 

The program presented in this paper is based on the methodics [8] where mixing coefficients ai L5 ^ (r) are 
transformed to a^^ (r) using the following relation: 

« s fo) (r) = Yl hsJPhr L r S r JP) (r), (5) 

r 

where (7 S JP\j r L r S r JP) is the LS — jj transformation matrix. 

In general, the performing of the transformations (5) with the matrix (j s JP\^ r L r S r JP) requires the 
LS — jj transformation matrices for the shell of equivalent electrons [10]. The numerical values of the 
matrices for the shell of equivalent electrons arc presented in [9] , nevertheless it is necessary to ensure the 
compatibility of the coefficients of fractional parentage (CFP) used to form (3) and (4) with the ones used in 
the transformation matrices. In this program the usage of the CFP obtained from the reduced coefficients 
(subcoefficients) of fractional parentage (RCFP) is assumed [1] (such CFPs were used for the calculation of 
LS — jj matrix elements). The program is based on quasisipin formalism and therefore the appropriate data 
are required. The data could be generated using the popular package [19, 20] with the angular parts [21]. 

3 Description of the Program 

The program is written in Fortran 90 programming language. The main part of the program consists of 
several subroutines collected into the module rabs_lsj j (file rabs_lsjj .f90). The program is designed in 
the way similar to Ratip [22] package. Some of the Ratip features are used and the overall program includes 
the additional modules rabs_angular, rabs_constant, rabs_csl, rabs_determinant, rabs_dirac_orbital, 
rabs_naglib, rabs_nucleus, rabs_rcfp, rabs_utilities, rabs_lsj, rabs_lsj_data [23, 24] (the corre- 
sponding files named as the module_name . f 90). 

3.1 Installation of the program 

The program was designed and tested on the PC with the Linux OS. Nevertheless it can be easily adopted 
to any Unix-Linux based platforms (including IBM RS/6000, Sun OS). The installation of the program is 
similar to Ratip components. First, the script make-environment which declares the global variables for 
the compilation and linkage of the program should be run via statement source make-environment. Then, 
the make-file (make-Is jj in this case) should be executed via comand make -f make-Is jj (both files can 
be obtained from the authors together with other program files). Then the executable xlsj j is generated. 

3.2 Interactive control and output of the program 

The principled scheme of the execution of the program is presented in Figure 1. After the initialization of 
several data arrays at the beginning of the execution (subroutine lsj_initialization_LS_j j () [23]), the 
control is taken by the procedure lsj j_control_transf ormationO , which reads and interprets all input 
data, carries out the transformation for all selected levels, and finally prints the results to screen. The 
example of an interactive dialogue performed through the execution is presented in Figure 2. 

3.2.1 Input 

The input of data in this version of the program is performed by the subroutine lsj j_load_csf s_f rom_MCHF_f iles () 
called from the subroutine lsj j_control_transf ormationO . The input data is assumed to be in the same 
format as the Mchf [19, 20] output data and can be generated with this package with the angular parts [21] 
(to ensure the usage of the CFP defined via the quasispin formalism). Two files, one with the CSF basis 
(configuration list file . inp) and the other with the ASFs expansions (mixing coefficients file . j ) should be 
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Table 1 The subroutines and functions of the rabs_lsj j module. 



General purpose 

lsj j_control_transf ormationO 
lsj j_load_csf s_f rom_MCHF_f iles () 
lsj j_form_csf _basis_j j () 
Is j j _transf ormat ion_ASF ( ) 
lsj j_get_subshell_term_j j () 

Output facilities 

lsj j_print_MCHF_data() 

lsj j_print_csf s_LS_MCHF() 

lsj j_print_coeff icients_LS_MCHF() 

lsj j_print_single_conf iguration_j j2() 

lsj j_print_coef f icients_j j () 

lsj j_print_conf iguration_scheme_j j () 

lsj j _print_detailed_inf o_csf _set_LS () 

Utilities 

lsj j_lval() 

lsj j_lvalr () 

lsj j_largest_coeff () 

lsj j_deallocate_asf _basis_j j () 



provided. The configuration list file may be generated using the Genclf and the mixing coefficients file - 
using the Cl program from the Mchf package [19]. The fragments of the input files are presented in Figures 
3, 4. The detailed specification of the format of these files is avialable at the description of the popular 
Mchf package (mentioned above) therefore we will not go into details. 

3.2.2 Output 

The output of the program, namely the jj-coupled CSFs and the expansions of the selected ASFs are 
performed either directly from the subroutine lsj j_control_transf ormationO or from the specific data 
output subroutines called from it. 

The direct "onscreen" output includes the quantum numbers and weights of several CSFs with the 
greatest weights as well as the sum of the squares of the mixing coefficients (in the output of the program 
denoted as Total sum of weights), serving as the numerical criteria of the accuracy of the expansions 1, 
2. 

The procedures lsj j_print_conf iguration_scheme_j j (), lsj j_print_coef f icients_j j () also print a 
full expansions of all ASFs in jj-coupled basis to the files lsj j-out-csf s . j j and lsj j-out-coef f . j j in 
the manner similar to the LS'-coupled data. The fragments of the output files are presented in Figures 5,6. 

3.3 Description of the module rabs_lsj j 

In this section we will introduce the main procedures of the module and explain some details of their usage 
because the main module of the program rabs_lsj j is assumed to be used within the Lsjj as well as in other 
programs (for example in other programs of the coupling optimization package Laris). The module collects 
the procedures and functions in an alphabetic order; they are briefly explained in the header of the module. 
Further information about the methods, data structures, etc. and their implementation can be found in many 
in-line comments in the headers of the individual procedures/functions or directly inside the source code. 
The names of the subroutines/functions in the module begin with the "lsj j_" and are followed by the rest 
of the names explaining the purpose of the procedure. Table 1 lists the procedures and functions grouped 
according to their purpose. The detailed description of the subroutines/functions is presented below. 
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The subroutine lsj j_control_transf ormationO 

The main subroutine is controling the flow of the execution. The subroutine performs an interactive dialogue 
with the user, reads the specification parameters for the calculation, opens and closes input and output 
streams, calls other subroutines of data input, the proceeding and the output. No direct input arguments 
should be specified. All the dinamically allocated arrays are deallocated at the end of the execution of the 
subroutine. 

The subroutine lsj j_load_csf s Jrom_MCHF_f iles () 

This subroutine loads the data on ASFs and corresponding expansions in terms of iS'-coupled CSFs from 
the files. The input files are assumed to be formated like MCHF output files [20]. Two input arguments 
stream_nr and stream_nr_j specify the input streams for the configuration list file and for mixing coefficient 
file. The output argument ierr is designed for error handling and is assigned to in case of failure. The 
data read are placed to the public variable asf _set_LS°/ csf _set_LS. 

The subroutine lsj j_f orm_csf _basis_j j () 

This subroutine generates the jj-coupled CSFs corresponding to the iS'-coupled ones stored in the public 
variable asf _set_LS°/ csf _set_LS. No direct input arguments should be specified. The data generated are 
placed to the public variable asf _set%csf _set. The algorithm of the generation is designed in a way that the 
subroutine can deal with the wave functions of the configurations with the arbitrary number of open shells. 
The algorithm is implemented using the utility of recursive subroutines. Two internal recursive subroutines 
are used. The first one lsj j Jorm_csf _basis_j j_f ill_occupations() is meant for the definition of the 

- + N 
number of jj-coupled shells (and occupation numbers Ni,Ni) for the given LS coupled shell l"\ And 

the second one lsj j_f orm_csf _basis_j j_job_count () is for the the calculation of the number of CSFs in 
jj-coupling and corresponding intermediate quantum numbers J- •, Ji..j (4). 

The subroutine lsj j_transf ormation_ASF() 

The subroutine transforms the atomic state function, represented in a LS-coupling CSF basis into a ba- 
sis of jj-coupling CSF. Namely the expansion coefficients a^^ (r) are calculated and placed to the array 
asf _set°/ asf (level_LS) "/eigenvector. The input argument level_LS specifies a serial number of the state 
to be considered. The LS — jj matrix elements ("f s JP\j r LS JP) for the transformation (5) are taken from 
module rabs_lsj_data via the subroutine lsj_transf ormation_LS_j j_general () [23]. 

The subroutine lsj j_get_subshell_term_j j () 

This procedure returns all allowed subshell terms (j, w, Q, J) for given j N which must be 1/2, 3/2, 5/2, 7/2 
or 9/2. Two input arguments namely j_shell and N define the value of the angular momentum and the 
occupation. The number of corresponding terms as well as their quantum numbers are given by the output 
arguments number and j j correspondingly. 

There are a number of subroutines for the data output in various forms. All the names for all subroutines 
begin with the lsj j_print_ and are followed by the rest of the names explaining the content and type of 
the output. 

The subroutine lsj j_print _MCHF_data() 

The subroutine prints the data in the manner of MCHF program. The input arguments streamnr_inp and 
streamnr_j specify the output streams for the MCHF mixing coefficients (. j) and configuration list (.inp) 
files correspondingly. The argument asf _set_LS specifies the variable where the LS'-coupled data are stored. 
The subroutine first generates the temporary set of CSFs not coupled to J (as it is in MCHF files), prescribes 
the mixing coefficients to this set of CSFs and stores this data in temporary variable asf _set_LS_MCHF. Then 
the subroutine calls lsj j_print_csf s_LS_MCHF() and lsj j_print_coef f icients_LS_MCHF() for the output 
of CSFs and ASFs with the corresponding mixing coefficients from this temporary variable. 

The subroutine lsj j_print_csf s_LS_MCHF() 
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The subroutine prints the configuration state functions to the file in the form similar to MCHF configu- 
ration list file (. inp). The argument streamnr specifies the output stream and the csf _set_LS specifies the 
CSFs in LS-coupling (not coupled to J). 

The subroutine lsj j_print_coef f icients_LS_MCHF() 

The subroutine prints the data on ASFs with expansion coefficients to the file in the form similar to 
MCHF mixing coefficients . j file. The argument streamnr specifics the output stream and the asf _set_LS 
specifies the ASFs and the mixing coefficients in LS'-coupling (for the list of CSFs not coupled to J). 

The subroutine lsj j_print_single_conf iguration_j j2() 

The subroutine prints all information about a single CSF in jj-coupling. The input argument stream 
specifies the output stream, the csf _set specifies the variable of type csf _basis where information on jj- 
coupled CSFs is stored and csf _number the serial number of CSF to be printed. The subroutine is similar 
to the lsj j_print_single_conf iguration_j j () defined in module rabs_lsj. 

The subroutine lsj j_print_coef f icients_j j () 

The subroutine prints expansion coefficients in j j-coupling to the file in the form similar to MCHF mixing 
coefficients (. j) file. The input argument streamnr specifies the output stream and the asf_set specifies 
the variable where the jj-coupled data are stored. 

The subroutine lsj j_print_conf iguration_scheme_j j () 

The subroutine prints information about the CSFs in jj-coupling. The input argument stream specifies 
the output stream and the csf _set specifies variable where the data on jj-coupled CSFs are stored. 

The subroutine lsj j_print_detailed_inf o_csf _set_LS() 

The subroutine prints detailed info about the LS-coupled CSFs. The input argument stream specifies 
the output stream and the csf _set_LS specifies variable where the data on LS'-coupled CSFs are stored. 

The module rabs_lsj j contains a number of auxiliary subroutines and functions. 

The function lsjj_lval() 

The function returns the integer value of an orbital quantum number given as a character. The range of 
an argument 1 should be from s to q (or from S to Q), i.e. the corresponding angular momentum I from 
to 12. 

The function lsjj_lvalr() 

The function returns the character value of an orbital quantum number given as an integer. The range 
of an argument 1 should be from 1 to 12. The subroutine is inverse to the lsj j_lval(). 

The function lsj j_largest_coef f () 

The function defines the serial number of the largest element in an array. The argument array specifics 
the dinamically allocated array and size - it's size. 

The subroutine lsj j_deallocate_asf _basis_j j () 

The subroutine deallocates the dinamically allocated data arrays within the variable asf _set (i.e. the 
jj-coupled data). The variable is defined in the module rabs_lsj. No direct input arguments should be 
specified. 

4 Example 

Here we present the example of the transformation of the atomic state function and classification in terms 
of quantum numbers of LS- and jj-coupling schemes. The example considers the excited state of Ne II. The 
interactive dialogue performed in the process of the execution of the program is presented in Figure 7. 

First the names of configuration list file (Ne . inp) and mixing coefficients file Ne . j ) are specified (the 
fragments of the input files presented in Figures 3 and 4). Then the program proceeds the input and asks 
to specify the ASFs to consider. 

After the specification of the ASFs the corresponding jj-coupled CSF basis is formed (144 CSFs in 
this case) and the mixing coefficients a s fo) (r) of the MCHF expansion (2) are calculated. The resul- 
tant CSF basis and mixing coefficients of the specified ASF are printed to the files lsj j-out-csf s . j j 
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and lsj j-out-coeff . j j correspondingly. The fragments of the output files are presented in Figures 5 
and 6. In addition the LS- and jj-coupled CSFs with the greatest weights (squares of the mixing co- 
efficients) useful for the identification and classification of ASFs under the consideration presented inter- 
actively. In our case the considered ASF may be identified by means of the quantum numbers of the 
LS-coupled CSF (((2s 2 ) 1 S, (2p 4 ) 3 P) 3 P, (Ss^S) 4 Ps\ or by means of the quantum numbers of jj-coupled 



CSF |((((2 S 2 )0,(2p- 2 )0) 0,(2p 2 )2) 2,(3s 1 )i) §>. 

The sums of the squares of the mixing coefficients in LS- and jj-couplings equal to 1 with the precision 
of the input coefficients in iS'-coupling (0.0000001). The fact that the sums in LS- and jj-couplings match 
up to the thirteenth digit reveals the accuracy of the performation of the transformation (5). 



5 Conclusions and outlook 

The program Lsjj enables us to transform the multiconfigurational expansions of the atomic state functions 
from LS- to jj-coupled bases. The input data in the current version of the program are supposed to be in the 
form of MCHF [20] output files, but the structure of the program enables to add more data input interfaces 
easily. The module rabs_lsjj together with the earlier constructed rabs_lsj [23] extends the possibility 
of comparison of ASFs expansions in various bases. This extent combined together with the program for 
the changing coupling scheme between the LS'-coupled shells of equivalent electrons and evaluation of the 
suitability of the intermediate quantum numbers for the identification and classification of atomic spectra 
under consideration [13] will serve as a tool for the search of the optimal set of quantum numbers for the 
classification of spectra of atoms and ions. 



Program is obtainable from Vilnius University Research Institute of Theoretical Physic and Astronomy, 
A. Gostauto 12, Vilnius, 2600, LITHUANIA. E-mail: tomas@mserv.itpa.lt. 
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PROGRAMA ATOMy BUSENy FUNKCIJOMS TRANSFORMUOTI IS LS I jj RYSIO BANGINiy 
FUNKCIjy BAZE 

T. Zalandauskas ir G. Gaigalas 
Santrauka 

Eksperimentiskai gaunami atomy ir jony energetiniai spektrai paprastai idcntifikuojami ir klasifikuojami 
naudojant LS rysio schemos kvantinius skaicius. Si rysio schema yra labai populiari ir teorineje atomo 
fizikoje, kai atomai ir jonai yra nagrinejami nereliatyvistiniame artutinumc. Bet LS rysio kvantiniai skaiciai 
daznai nera tinkami sunkiy atomy bei didelio ionizacijos laipsnio jony busenoms apibudinti. Vienareiksmiam 
bei tiksliam tokiy atomy bei jony buseny idcntifikavimui bei spektry klasifikavikavimui geriau tinka kity rysio 
shemy (daznai jj rysio) kvantiniai skaiciai. 

Neseniai gautos LS-jj transformacijos matricos ekvivalentiniy elektrony sluoksniams (iskaitant / elek- 
trony sluoksnius su bet kokiu uzpildymo skaiciumi) bei LS-jj transformacijos matricy israiskos konfiguracijoms 
su bet kokiu atviry sluoksniy skaiciumi igalino mus sukurti programa, skirta atomy bei jony buseny funkcijy 
transformavimui is LS i jj rysio banginiy funkcijy baz§. Tokio pobudzio programy, leidzianciy keisti rysio 
tipa tiek ekvivalentiniy elektrony sluoksniy viduje (iskaitant / elektrony sluoksnius), tick ir tarp sluoksniy, 
iki siol nebuvo. 

Programa igalina atlikti atomy ir jony spektry identifikavima ir klasifikacija LS ir jj rysio kvantiniais 
skaiciais, praplecia kvantiniu skaiciu rinkinio, optimalaus atominiu spektry idcntifikavimui ir klasifikavimui, 
paieskos galimybes. 
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I lsj_initialization_LS_j j I 
| 1 



I lsj j_control_transf ormationO I 
| 1 



data input 

- lsj j_load_csf s_from_MCHF_f ilesQ 



formation of all CSFs in jj coupling 

- lsj j_f orm_csf _basis_j j () 



LS-jj transformation and classification 

of all selected levels 

- lsj j_transf ormation_ASF() 



data output 

- lsj j_print_conf iguration_scheme_j j () 

- others 



Figure 1 : Schematic flow chart of the Lsjj program. 

Program for transformation of atomic states from LS to jj coupling 
(C) Copyright by T Zalandauskas and G Gaigalas Vilnius (2003) . 
Transform one or several ASF from a MCHF calculation into a 
jj-coupled CSF basis. The transformation starts from the given 
cfg.inp and .j files and is carried out for the n leading CSF in 
the LS-coupled basis; the new representation in the jj basis is 
printed as in a standard GRASP92 computation. 

Enter the name of the MCHF configuration list file: 

Enter the name of corresponding .j mixing coefficient file: 

load_cfg_from_MCHF_f ile . . . 
number of core shells = ... 
total number of shells = ... 
number of csf s LS = . . . 

reordered 

there are . . . nonrelativistic CSFs (before couplin to J) ; 
now couple them to J ... 

now load mixing coefficients and form asf _set_LS . . . 
end reading j 

there are . . . atomic state functions 
there are . . . nonrelativistic CSFs (coupled to J) 
. . . load complete . 
Maximum number of considered ASF is : ... 

Enter the level numbers of the ASF which are to be transformed, 
e.g. 1 3 4 7 - 20 48 69 - 85 : 



P 
R 

G 
R 
A 
M 



E 
X 
E 
C 
U 
T 
I 

N 



Figure 2: Interactive dialogue of Lsjj program 
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Figure 3: Fragment of MCHF configuration list file Ne.inp 



Ne Z = 10.0 N = 9 NCFG = 224 



2*J = 5 NUMBER = 1 
3 -126.95806136 

0.0000000 0.0000000 0.9951298 0.0130359 0.0003521 0.0000000 0.0237319 

0.0096379 -.0179847 -.0159071 0.0000000 -.0002664 0.0000000 0.0177844 

0.0000002 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 

0.0575015 -.0000647 -.0005302 0.0000000 0.0001279 0.0000000 0.0000000 



Figure 4: Fragment of MCHF mixing coefficients file Ne. j 
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The current configuration scheme with 144 CSF in jj coupling 

is defined as follows: 



Number of electrons : 9 

Number of (relativistic) subshells: 9 
Core shells : 
Is 

Peel shells: 

2s 2p- 2p 3s 3p- 3p 3d- 3d 

144 CSF(s): 

1) 2s ( 2) 2p-( 1) 2p (3) 3s ( 1) 

1/2 3/2 1/2 

1/2 2 5/2 

2) 2s ( 2) 2p-( 2) 2p (2) 3s ( 1) 

2 1/2 

2 5/2 

3) 2s (2) 2p ( 4) 3d ( 1) 

5/2 

5/2 

4) 2s ( 2) 2p-( 1) 2p (3) 3d ( 1) 

1/2 3/2 5/2 

1/2 1 5/2 

5) 2s ( 2) 2p-( 1) 2p (3) 3d ( 1) 

1/2 3/2 5/2 

1/2 2 5/2 

6) 2s ( 2) 2p-( 1) 2p ( 3) 3d-( 1) 

1/2 3/2 3/2 

1/2 1 5/2 



Figure 5: Fragment of output file lsj j-out-csf s. j j of jj coupled CSFs 



1 -126.95806136 

0.5638947 0.8200464 -.0140585 0.0254398 0.0123155 -.0067862 -.0061464 

0.0103721 0.0170042 -.0084857 0.0230858 -.0169735 0.0440080 -.0234300 

0.0355732 0.0514159 -.0002857 -.0018696 0.0043991 0.0004125 -.0016268 

0.0049082 0.0044622 -.0028572 -.0023440 0.0025700 0.0059582 0.0018875 



Figure 6: Fragment of output file lsj j-out-coef f . j j of the ASFs expansions in jj-coupling 
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Enter the name of the MCHF configuration list file: 
Ne . inp 

Enter the name of corresponding .j mixing coefficient file: 
Ne.j 

load_cfg_from_MCHF_f ile . . . 
number of core shells = 1 
total number of shells = 6 
number of csfs LS = 224 

there are 224 nonrelativistic CSFs (before coupling to J) ; 
now couple them to J ... 

now load mixing coefficients and form asf _set_LS . . . 
end reading j 

there are 1 atomic state functions 
there are 144 nonrelativistic CSFs (coupled to J) 
... load complete. 
Maximum number of considered ASF is : 1 

Enter the level numbers of the ASF which are to be transformed, 
e.g. 1 3 4 7 - 20 48 69 - 85 : 

1 

number_of _levels , levels(:) =11 
start formation of CSFs in jj coupling. . . 
number of shells in jj coupling = 9 
number of core shells in jj coupling = 1 

number of configuration state functions in jj coupling = 144 
finish formation of CSFs in j j coupling. . . 

Weights of major contributors to ASF in LS-coupling: 
Level J Parity CSF contributions 

1 5/2 + 0.99028 of 1 0.00391 of 15 0.00331 of 11 

Definition of leading CSF: 

1) 2s( 2) 2p( 4) 3s( 1) 

ISO 3P2 2S1 3P 4P 5/2 

Total sum of weights is: 0.9999999614146498 



Weights of major contributors to ASF in jj -coupling: 
Level J Parity CSF contributions 

1 5/2 + 0.67248 of 2 0.31798 of 1 0.00264 of 16 

Total sum of weights is: 0.9999999614146532 

Definition of leading CSF: 

1) 2s ( 2) 2p-( 1) 2p (3) 3s ( 1) 

1/2 3/2 1/2 

1/2 2 5/2 

2) 2s ( 2) 2p-( 2) 2p (2) 3s ( 1) 

2 1/2 

2 5/2 
16) 2s ( 1) 2p-( 2) 2p (2) 3s ( 2) 
1/2 2 

1/2 5/2 5/2 

LS-j j complete ... . 



Figure 7: Interactive dialogue of LSJJ program for the transformation of state function of excited Ne 
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